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We consider spin ice magnets (primarily, Dy2Ti2 0?) in the vicinity of their critical point on the (H,T) 
plane. We find that the longitudinal susceptibility diverges at the critical point, leading to the behaviour 
qualitatively similar to the one which would result from non-zero conductance of magnetic charges. We show 
that dynamics of critical fluctuations belongs to the universality class of easy-axis ferroelectric and calculate 
logarithmic corrections (within two-loop approximation) to the mean-field critical behavior. 

PACS: 75.40.Gb, 75.47.Lx 



Spin ice [T] is a kind of geometrically frustrated mag- 
netic materials with a structure of low-energy states 
presenting deep similarity with usual water ice, those 
H-bond structure was analized long ago by Pauling [2]. 
The most studied examples of spin ice include py- 
rochlorc oxides Dy2Ti2C>7 and Ho2Ti207; their struc- 
ture is shown in Fig.Q] Spin ice demonstrates a number 
of peculiar low-temperature properties, the major of the 
are: i) an extensive residual entropy, and ii) elementary 
excitations resembling magnetic monopoles [21 [T]. Low 
temperature dynamics of spin ice is governed almost 
solely by these magnetic excitations, therefore one may 
expect spin ice to demonstrate phenomena similar to 
those known for the electrolytes. Indeed, it was recently 
shown that some aspects of spin ice behaviour can be 
described in terms of "magnetolyte" similar to the On- 
sager theory [3] of electolytes [5]. It is also known pQ, 
that at low enough temperatures spin ice undergoes a 
first-order transition as function of applied magnetic 
.field H . The first-order transition line in the (H, T) 
plane terminates at the critical point H Cl T c . In particu- 
lar, H c = 0.9T and T c = 0.38K in the case of Dy 2 Ti 2 7 . 
In the present Letter we calculate dynamic longitudi- 
nal magnetic susceptibility of spin ice near this critical 
point; our results predict the behaviour formally simi- 
lar to the one expected for the media with a nonzero 
"magnetoconductance" . We will explain, however, that 
'it would be incorrect to interpret these results in terms 
of "direct current" of monopoles. Specific numerical es- 
timates will be done for Dy2Ti2C>7, whereas our general 
scheme is applicable to any pyrochlore spin ice. 

The crucial difference between spin ice and elec- 
trolyte is the finite static magnetic susceptibility of 
the former as opposed to the diverging dielectric low- 
frequency response e(w) cx icr/co of any electrolyte with 
a nonzero conductivity er. The origin of this difference 
can traced to the fact that magnetic monopoles of spin 




Fig.l. Pyrochlore lattice consists of corner-sharing 
tetrahedra; spins of dysprosium are siting at their ver- 
tices. 

ice are the sources of magnetization field M (or, equiva- 
lently, magnetic field H), but they do not generate mag- 
netic induction B, which is the proper analog of electric 
field E for electrostatics [5]. 

Following [7] , one can represent linear magnetic sus- 
ceptibility of spin ice as 



Xo\T)-iw/T 



(1) 



where Xo(T) ~ ^l 2 b /{VqT) is static susceptibility ( /is is 
the Bohr magneton and Vq is the volume of an elemen- 
tary cell). Therefore relaxation of magnetization after 
the sudden switching of magnetic field is given by 



M(t) = X oH(l - exp[- 



(2) 



At the small timescales t <C Xo/r the relaxation pro- 
cess (0) resembles direct conductivity with "magnetic 
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Fig.2. Pyrochlore lattice can be represented as a stack 
of triangular and kagome lattices, with easy axis di- 
rections being aligned with the bonds of diamond lat- 
tice, which is dual to the pyrochlore lattice. Magnetic 
monopoles are situated in the centers of tetrahedra. 

current" proportional to the external magnetic field: 
j m = dM/dt = TH, thus the magnetic analog of con- 
ductivity is cr jn = r. However, finite value of static 
susceptibility Xo(T) limits possible observation of mag- 
netoconductivity to a transient regime only. 

We emphasize that Eq.((T]) for dynamic susceptibil- 
ity xo(H, T) with some effective relaxation rate T(H, T) 
should be valid in a broad range of temperatures and 
fields, when monopoles are present. Spin ice in the ab- 
sence of external field, and spin ice in sufficiently strong 
field when kagome ice structure develops, represent the 
two well-defined limiting cases. Below we consider the 
vicinity of the critical point [H C ,T C ) where static sus- 
ceptibility is expected to diverge, leading to a broad 
timescale for the observation of magnetoconductivity. 
We will use virial expansion in order to estimate the 
parameters of the Ginzburg-Landau free energy F{m} 
which describes fluctuations in the critical region. Then 
we use this F{m} functional to construct an effective 
Martin-Siggia-Rose functional [8] HI [10] which describes 
critical dynamics of spin ice, and study it up to the two- 
loop logarithmic approximation. 

The Hamiltonian of the spin ice is of the form[T] 

H = -~y] Jjjjej ■ ej)SiSj - ^sH^ejSj (3) 



<\ • <\ 3(e, • Tij)(ej ■ Tij) 



Si Sj 



Here Si = ±1 is an effective Ising variable which de- 
scribes sign of magnetic moment fii = S^ei of the i- 
th Dysprosium ion; the magnitudes \fi{\ = fi « 10/is, 
whereas are the unit vectors along easy axis, see 



Fig. [21 the direction of axes being chosen with respect 
to one of the tetrahedral sublattices, see Fig.[TJ Param- 
eters Jij describe exchange interactions between first-, 
second- and third-order neighbours, their magnitudes 
were determined [11] as J\ = — 3.72K, J 2 = 0.1K and 
J3 = — 0.03K respectively. External magnetic field H is 
directed along the [111] axis, and the last term in Eq.Q 
stands for the magnetic dipole-dipole interaction, the 
dipole interaction constant D = /j, 2 /a 3 = 1.41K, where 
a = 3. 54 A is the pyrochlore lattice constant (it is equal 
to the nearest-neighbour distance, see Fig. [2]). 

We will study the system described by the Hamil- 
tonian ([3]) near the critical point, where large average 
magnetization My along the field direction is induced. 
Deviation m of the actual magnetization My from its 
average value at the critical point M^(H C ,T C ) can be 
considered as an order parameter which describes the 
state of the system in the vicinity of the critical point. 
The corresponding Ginzburg-Landau free energy func- 
tional is of the form 



dV 
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(4) 



where a = a(T—T c ) and h = H— 7(T— T c ). In a system 
with local interactions one would find b = bk 2 = — 6A, 
but our case is different due to the presence of strong 
dipole-dipole interactions. Similar problem of second- 
order phase transition in uniaxial ferroelectric was con- 
sidered in the seminal paper (12) . where it was shown 
that in presence of dipole-dipole interaction operator b 
should be modified as follows: 



b = bk 2 + 4ttx 



(5) 



where x = k z /k and the z is [111] axis. We are inter- 
ested here in the frequency-dependent susceptibility of 
spin ice, and thus we need to extend the renormaliza- 
tion group analysis developed in Ref. [12] for the criti- 
cal dynamics. Due to absence of any locally conserved 
quantities in the problem, the system can be described 
by purely relaxational dynamics governed by the free 
energy (g]): 

an.~r£ (0) 

To obtain estimates for coefficients in Q we use virial 
expansion as described below; the estimate for the ki- 
netic coefficient V will be provided at the end of this 
Letter. At low temperatures T <§C T c , in the high-i7 
phase almost all spins are aligned with the magnetic 
field, while in the \ow-H phase the same is valid for 
the spins of triangular sublattices only; configurations 
of other spins is governed by the kagome ice rules, see 
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Fig. 2 and Ref. [T]. One possible way to estimate the 
parameters of T near the critical point would be to con- 
sider gas of interacting monopolcs taking into account 
their direct as well as entropic interactions |13) : similar 
problem of lattice ion systems was considered in [14]. 
However, we prefer to start from the high-i? phase and 
to consider an exponentially dilute (far from the critical 
point) gas of flipped spins as a lattice gas of interacting 
particles. Such an approach will be approved a poste- 
riori by the numerically small concentration of flipped 
spins even at the critical point. Using the Hamiltonian 
([3]) we obtain the energy of such a particle (siting in the 
site 0) in the form 



E = -/iff 



2Y^ Joi(eo • ei)Si - 



(7) 



e -e.; 3(e • r 0i )(ei • r 0i 



|r ; 
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where S% = 1 on triangular sublattices and Si - 
on kagome sublattices. The interaction of two particles 
siting in the sites i and j is equal to 



Uij = -Jij +4Da 3 



(8) 

Now we can employ virial expantion |15j in terms of 
small density of particles n, to obtain free energy 

E 



F(n) = TN [~l-n-b 2 n 2 -{b 3 -2bi)n J 



T 



- n In n] 
(9) 



U. 



where T-dependcnt functions 

h = ^ ^2 foi, with fij = cxp 

i 

^ 3 = 31 ^[/ 0i Aj'^0 + hifoj + fiofij + foifoj] 
i J 

are "cluster integrals" given by the summation over 
the lattice. We calculate these sums numerically and 
find the position of the critical point H^ h = 1.37T, 
T*" = 0.85-ftT and critical concentration of "particles" 
n c = 0.14 from the conditions 

dF d 2 F d 3 F , . 

^T = °> TT = ^' TT = ° ( 10 ) 
dn dn z dn^ 

Smallness of n c supports qualitative validity of our 
virial expansion; however, experimental values of H c 
and T c are below our estimates by factors 1.5-2. With 
the position of the critical point being determined, we 
can estimate values of the parameters entering free en- 
ergy functional T: 

da 



dT 
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Fig. 3. Diagrams responsible for the renormalization 
of coefficients 6 and F Solid lines represent correla- 
tion function (mim2), and mixed lines correspond to 
response function (mip2). 

However, virial expansion is not useful to determine the 
coefficient b entering the gradient term in Q; thus we 
estimate it within nearest-neighbours approximation, 
neglecting long-distance tail of the dipole-dipole inter- 
action: 

5V2Ji + 5L> 2 _ S 2 



26A' 



(12) 



24 D 

Coming back to the problem of longitudinal suscep- 
tibility and considering it in the mean-field approxima- 
tion, one obtains (1) with Xq 1 — a (T — T c ) : so exactly 
in the critical point \ ^ ir/cj. Wc will sec below that 
account of fluctuations do not change this result consid- 
erably, leading to logarithmic corrections only. 

With the parameters of the free energy @ speci- 
fied, we move to the analysis of critical fluctuations and 
their role in the dynamics. Thus wc construct, using 
Eqs. pio) . the Martin-Siggia-Rose [HI El HO] dynamic ac- 
tion 

X = [ dtdV\r~ 1 pdtm+apm+pbm-\ — prr? — ph+iTT^ 1 p 2 ] 
J 6 

(13) 

where p = p(t, r) is the dynamic field conjugated to the 
magnetization field m — m(t, r), and the last term in 
the action describes thermal noise. The estimate for the 
kinetic coefficient T will be provided later. Summation 
of leading logarithmic corrections with the action (TIB"]) 
up to one-loop approximation can be done in complete 
analogy with the paper of Larkin and Khmelnitsky [12j , 
as it refers to the static quantities only: 



dX da 1 



ST 
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?= 2 ln Z^ 



(14) 



(15) 



where A 



is an ultraviolet cutoff. Renormaliza- 



tion of the gradient term b and the kinetic coefficient 
r appears in the second-loop approximation only. The 
corresponding diagrams are represented in FigJ3] Cal- 
culation of these diagrams leads to the renormalization- 
group equations 
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where the constant C 
lfj"6l) is as follows: 



400. The solution of equations 
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r = Tq exp 



In- 
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3 3 7T 6 

6A 2 



bk 2 + 4irx 2 



dpds direction (1,1,1) of the applied field. Thus positive- 

charge and negative-charge monopolcs will be siting, 
preferably, on the "upper" and "lower" sublattices cor- 
respondingly, with the energy difference /iH/3 between 
them. As a result, magnetization is primarily deter- 
mined by the concentration of the monopoles. 

Finally, we need to estimate relaxation rate T. Dy- 
namical processes in spin ice at low temperatures and in 
low-field phase arc governed almost solely by monopoles. 
Ground state doublet of each spin is separated by large 
energy gap A ~ 300K from higher-energy states([T|). 
Thus, at temperatures T < IK, processes of spins flips 
are solely quantum. Then, the dependence of magneto- 
conductivity on temperature follows Arrhcnius law with 
activation energy E m ss IK corresponding to the cre- 
ation of one monopole. Reminding the relation between 
r and <7, we find 



(17) 



Using Eqs. (|llll5p we find that initially value of coupling 
constant is go = 0.15 <C 1. Then renormalization of b 
and r is very small numerically and we neglect it, so the 
susceptibility takes the form 



1 



X 



a(k, x, uj) + bk 2 + 47T.T 2 — ico/T 
and in the uniform field 



(18) 
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X = 



ico/T ' 



o(T) = 



a(T - T c ) 



1 



1 6A 2 

9o m a \ T -T c 



+M/r; 



(19) 

For T < T c static susceptibility should be modified by 
a factor of two: a should be substituted by a/2. This 
is the result in the absence of effective field h which is 
defined after Eq.Q. If the temperature T = T c , then 
the renormalization cutoff will be determined by h, i.e. 
by the deviation H — H c 



1 



iui/T ' 



a*(h) 



%\oh 2 Y 



(1 + ig In 



6A 2 



(A fi 2 )3+M/r 



(20) 

Exactly at the critical point a = and x ~ * w ■ 
However, in this system there is no conductivity in the 
usual sense. The point is that our lattice consists of a 
stack of interchanging kagome and triangular layers. In 
moderate magnetic fields (in particular, in the vicinity 
of the critical point), orientation of spins on the triangu- 
lar layers are almost fixed. This results in the situation 
when monopoles are bound to theirs kagome layers and 
thus there is no direct motion of monopoles in the direc- 
tion of the field. In fact, monopoles of one kagome layer 
reside in centers of tetrahedra, i. e. on two sublattices 
that are separated by a finite distance h = ad/5 in the 



T(T,H = 0)=Toe 



-E m /T 



(21) 



where Tq is temperature- independent. Low- field dy- 
namical magnetic susceptibility of Dy2Ti2C>7 was mea- 
sured in Ref. |TB] . Using the results [TB] for the imagi- 
nary part of susceptibility together with Eq.([T]), we ob- 
tain 

r(T = 1.8K, H = 0) ~ 100s" 1 

Comparing it with Eq. (|2ip and the estimate for E m . we 
conclude that To ~ 100s _1 as well. 

The nature of spin relaxation in spin ice near critical 
point is not quite clear. As we already mentioned above, 
in presence of magnetic field ~ IT, positive and negative 
magnetic charges reside on different sublattices and the 
tunneling processes become essentially modified, as the 
positions of monopolcs on the kagome lattices are not 
equivalent. There are two possible scenarios for tunnel- 
ing processes: sequential and simultaneous flips of two 
spins. Below we estimate the rate of sequential pro- 
cess, which will provide the lower bound for the relax- 
ation rate near critical point. Energy difference for the 
monopolcs of the same sign siting in two kagome sub- 
lattices leads to the exponential ratio of their concen- 
tration: N 2 = Nxexp[—2nH/3T\. On the other hand, 
detailed balance condition reads iViFi = N2T2 where 
subscript "1" refers to monopoles situated on "native" 
sublattice, and subscript "2" to "foreign" ones. Assum- 
ing T2 ~ To since no additional barrier exists for the 
hopping of monopole from a "foreign" site to a "native 
one" , we come to the estimate 

Y u PP er (H, T) - r (0, T) cxp[-2 i u J ff/3T] (23) 

and eventually T c - (0.01 ~ 100)s -1 . 
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In conclusions, we calculated dynamic spin suscep- 
tibility of spin ice material Dy2Ti2 07 near the critical 
point. Logarithmic corrections due to critical fluctua- 
tions are found to be small. The response of the form 
x(uj) ~ (iuj)^ 1 is predicted in the broad range of low 
frequencies. Measurement of the prefactor in this de- 
pendence would allow to determine the nature of ele- 
mentary spin flip processes near critical point. 
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